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Abstract 

Starting from a potential with a continuum of energy eigenstates, we show how the 
methods of supersymmetric quantum mechanics can be used to generate families of 
potentials with bound states in the continuum [BlCs] . We also find the corresponding 
wave functions. Our method preserves the spectrum of the original potential except 
it adds these discrete BlCs at selected energies. Specifically, we compute and graph 
potentials which have bound states in the continuum starting from a null potential 
representing a free particle and from both the attractive and the repulsive Coulomb 
potentials. 
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I. Introduction 



In 1929, Von Neumann and Wigner realized that it was possible to construct 
potentials which have quantum mechanical bound states embedded in the classical en- 
ergy continuum (BICs). Further developments, by many authors [2-6] have produced 
more examples and a better understanding of the kind of potential that can have 
such bound states, although there is not as yet a fully systematic approach. These 
authors have also suggested possible applications to atoms and molecules. Further- 
more, these works have shown that such BICs appear mainly in certain oscillatory 
potentials whose envelopes fall off fast enough to lead to normalizable wave functions 
but sufficiently slowly, such that the different maxima are able to conspire to create 
a captive state. Friedrich and Wintgen have given the example of two conspiring 
resonances and also of a hydrogen atom in a uniform magnetic field. Robnik |^ has 
shown in a similar way that a simple separable Hamiltonian can develop bound states 
in the continuum. In his examples too, coupled channels are responsible for the cre- 
ation of the BIG. Such a BIG is a very fragile structure - a small perturbation of the 
potential transforms it into a decaying resonance. Nevertheless, Gapasso et al. 
have recently reported direct evidence for BIGs by constructing suitable potentials 
using semiconductor heterostructures grown by molecular beam epitaxy. Finally, it is 
interesting to note that BIGs have found their way into a text IQ , illustrating for stu- 
dents the surprising possibility of the existence of quantum mechanical bound states 
in the classical continuum. 



Recently, extensive work has been devoted to generating isospectral potentials by 
the methods of supersymmetric quantum mechanics (SUSYQM) [9 - 14]. Starting 
from the Schrodinger equation for a potential, whose ground state wave function is 
known, this method permits one to generate families of new potentials, which may 
look quite different from the original one, but have exactly the same spectrum [12-14]. 
These methods are based on procedures invented by Darboux |jl5[ and generalized 
by Grum [T^. In this paper we extend the usual SUSYQM formalism for obtaining 
isospectral potentials and apply it to potentials with a continuum of scattering states 
to generate new potentials with bound states in the continuum. We show that, while 
the wave functions in the continuum of the original potential are non- normalizable, 
the ones generated by SUSYQM are normalizable thus representing a bound state. 
In particular, we construct one-parameter and two-parameter families of supersym- 
metric partner potentials with one and two bound states in the continuum. 



2 



In Section II, we present the development of the one-parameter family in SUSYQM 
to generate potentials with a single bound state in the continuum. The procedure is 
valid for any spherically symmetric potential V(r) which vanishes a.s r oo. We il- 
lustrate our procedure with two explicit examples: (A) free particle = on the half 
line; (B) attractive and repulsive Coulomb potentials. In Section III, we generalize 
the SUSYQM method to the two-parameter family of isospectral potentials and show 
that it generates two bound states in the continuum. Section IV contains a summary 
of our results. The method can readily be extended to an arbitrary number of BICs. 



II. The One Parameter Family of BICs. 

The radial s-wave Schrodinger equation for the reduced wavefunction u(r) (in units 
where h =2m=l) is 

-u" + V{r) u{r) = E u{r), (1) 

where we have scaled the energy and radial variables such that all quantities are 
dimensionless. A prime denotes differentiation with respect to r. For any potential 
which vanishes at infinity, Eq. (1) has a classical continuum of positive energy solu- 
tions which are clearly not normalizable. 

Using the formalism of SUSYQM and the Darboux |15| procedure for deleting and 



then reinstating the ground state Uo{r) of a potential V(r), one can generate a family 
of potentials V{r;X) which have the same eigenvalues as V(r). These isospectral 
potentials are labeled by a real parameter A which lies in the ranges A>OorA<— 1. 
The isospectral potential V{r] A) is given in terms of the original potential V(r) and 
the original ground state wave function Uo{r) by [12,14] 

V{r; A) = V{r) - 2[/n(/o + A)]" = V{r) - |^ + ^j^,, (2) 

where 

/o(r)^ rul{r')dr'. (3) 







Let us recall the main steps in obtaining Eq. (H). First, one writes the Hamiltonian 



H = —-^ + V{r) in factorized form H = A^A, with the operators 



A=^ + W{r), A^ = -^ + W{r) 
dr dr 
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The superpotential is given by W{r) = —Uq/uq. The supersymmetric partner Hamil- 
tonian is 



where 



u'r 



V+{r) = W^ + W' = V{r) -2[^ 



If the potential V(r) has eigenf unctions u„(r) at energies then the SUSY partner 
potential V+{r) has the same energy eigenvalues as V(r) with eigenfunctions AM„(r), 
except that there is no ground state at E=0 since Auo{r) = 0. This is the standard 
procedure for deleting the ground state and obtaining V"+(r). To re-insert the ground 
state, one asks for the most general superpotential W{r) such that 

VJr) = W^ + W' 



and this can be shown to be 17 



W{r; A) = W{r) + ^In {Io{r) + A), 
ar 

with Iq given in Eq. @. Thus the entire family of potentials V{r; X) = W^lr; A) — 
W'{r; A) has the same supersymmetric partner potential V+(r) obtained by deleting 
the ground state. 

In all previous work, uq was taken to be the nodeless, normalizable ground state 
wave function of the starting potential V(r). However, for the purposes of this paper, 
we can generalize the above equations to the case where uo{r) is any solution of 
Eq. (1) with arbitrary energy Eq. If Uq^v) has nodes, this leads to singular super- 
potentials and to singularities in the partner potential V+{r). However, when the 
original state at Eq is re-inserted, the resulting family of potentials V{r; A) is free of 
singularities [0. Our results are best summarized in the following statement: 

Theorem: Let Uo{r) and Ui{r) be any two nonsingular solutions of the 
Schrodinger equation for the potential V(r) corresponding to arbitrarily selected en- 
ergies Eq and Ei respectively. Construct a new potential V{r] A) as prescribed by Eq. 
(0). Then, the two functions 

UQ{r;X) = ^y (4) 
iQ -r A 

and 

Mi(r; A) = (£"1 - Eo)ui + Uq W{uq, Ui), (5) 
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[where W denotes the Wronskian, W{uo,ui) = uou[ — uiUq] are solutions of the 
Schrodinger equation for the new potential V{r; X) corresponding to the same ener- 
gies Eq and El. 

While the new potential in Eq. (|^) and the new wave functions in Eq. (^) were 
originally inspired by SUSYQM, the easiest proof of the above theorem is by direct 
substitution. One simply computes —it'- + V{r; \)ui (i=0,l), with the wave functions 
Ui given in the theorem. After straightforward but tedious algebraic manipulations, 
one gets Eiiii, thus establishing the theorem. The algebra is considerably simplified 
by using the following identity for the Wronskian of two solutions of the Schrodinger 
equation: 

— W{Uo, Ml) = {Eq - Ei)UoUi. (6) 

ar 

In the present work, we will take uq to be a scattering solution at a positive energy 
Eq = fc^ of a potential V(r) which vanishes at r=oo. Taking Mo(r = 0) = satisfies 
one of the required boundary conditions, but clearly uq oscillates as r — oo and has 
an amplitude which does not decrease. Consequently, the integral Io{r) in Eq. @ 
grows like r at large r and Uq is now square integrable for A > 0, while the original 
wave function uq was not. Therefore, we see that all the potentials V{r; A) have a 
BIG with energy Eq. Note from Eq. that -uq has the same zeros as the original 
Uq. At zeros of Uq, V{r; A) and V{r) are equal. On the other hand, all other oscilla- 
tory solutions to V{r) get transformed into oscillatory solutions to V{r; A) with the 
same energy. In particular, note that Ui{r; A) remains a non-normalizable scattering 
solution of the corresponding Schrodinger equation. 

We note that the new potential V{r; A) in Eq. (0) and the BIG at energy Eq are 
formed using the corresponding wave function UQ{r). Any other state, say Ui{r), is 
transformed into a solution of the new Schrodinger equation by the operation given 
in Eq. which involves both uq and Ui. The central column of the table gives a 
convenient overview of the relationship of the potentials V and V and the solutions 
of the corresponding Schrodinger equations. 

We now give two examples to explicitly illustrate how one applies the above pro- 
cedure to obtain potentials possessing one BIG. 

A. Free Particle on the Half Line. 
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Here, we consider the case V = 0, the free particle on the half line < r < oo. We 

choose Mo = sinkr, the spherical wave solution, corresponding to energy 

Eq = k'^ > 0, which vanishes at r = 0. The integral Iq given in Eq.(^ becomes 

Jo = [2kr - sin{2kr)]/{4:k). (7) 

We observe that /o^r/2asr^oo. 

The potential family V, defined in Eq.(0) becomes 

32 fc^ sin'^kr 8 k'^sin(2kr) , , 

''<'-'^) = ^5^ 

with 

Do(r; A) = 2kr - sin{2kr) + AkX. (9) 
V has a BIC at energy Eq = k"^ with wave function 

Mo(A) = 4A; sinkr /Dq. (10) 

For special values of the parameters k and A, the potential V and its BIC wave 
functions are shown in Figs, la and lb. The original null potential has now become 
an oscillatory potential which asymptotically has a 1/r envelope. The new wave func- 
tion at Eq = k"^ also has an additional damping factor of 1/r which makes it square 
integrable. As Uq appears in the numerator of V, Eq. (g), every node of Uq is asso- 
ciated with a node of V but not every node of V produces a node of Uq. The value 
of the eigenenergy Eq clearly is above the asymptotic value, zero, of the potential. 
Evidently, the many oscillations of this potential, none of them able to hold a bound 
state, conspire in such a way as to keep the particle trapped. 



The parameter A which appears in the denominator function Do{r; A) plays the 
role of a damping distance; its magnitude indicates the value of r at which the mono- 
tonically growing integral Iq becomes a significant damping factor, both for the new 
potential and for the new wave function. This is illustrated graphically in Figs, la 
and lb which are drawn for very different values of A. [Note that the wave functions 
shown in the figures are not normalized]. The parameter A must be restricted to 
values greater than zero in order to avoid infinities in V and in the wave functions. 
In the limit X ^ oo, V becomes identical to V. 
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B. Coulomb Potential 



Starting from the potential V = Z/r, for either positive or negative values of Z, one 
can easily construct the one-parameter family of isospectral potentials possessing a 
normalizable positive energy wave function. Here the unbound, reduced 1 = wave 
function satisfies the Schrodinger equation Eq.(l), which can be written in standard 
form 

4' + (1 - 2^)mo = (11) 

with p = \fEr and r] = Z/2\fE. 



For both positive and negative f/, the solutions involve confluent hypergeometric 
functions which in the asymptotic limit approach sine waves phase-shifted by a loga- 
rithmic term. Useful expressions for these solutions in the regions near and far from 



the origin are available in the literature |19, 20|. Stillinger and Herrick 0, following 
the method of Von Neumann and Wigner [|l| , have constructed BIC potentials and 
wave functions for the case of the repulsive Coulomb potential. Here we use our the- 
orem to construct a one-parameter family of isospectral potentials containing a BIC. 
The procedure is the same for both positive and negative Z; the only difference being 
in the sign of fj. The formal expressions for the BIC potentials and wave functions 
have been given above, Eqs.(0) and (0), in terms of uq. 



For both the attractive and the repulsive Coulomb potentials, the positive energy 
solution of Eq. (|l^) can be written in the usual form [|l^, |^ as the real function 

Uo{p) = Coiv) e-'m{l - If,, 2, 2tp), (12) 

where 

Co{f,) = {e-^'^'')\T{l + if,) I (13) 

and M(a,b,z) is Kummer's function. Using tabulated expressions for the Coulomb 
wave functions |^ and doing the integral for Jq numerically, we have obtained the BIC 
wave functions for representative values of A. The corresponding one-parameter fam- 
ily of potentials obtained by the SUSY procedure is given in Eq. (H) with Vo = Z/r. 



The results are displayed in Figs. 2 and 3. Fig. 2a shows the BIC partner to 
the attractive Coulomb potential for A = 1, k = 1, and Z = -2. Fig. 2b shows the 
(unnormalized) wave function of the bound state in the continuum for this potential 
at -Eq = k^- Fig- 3a shows the BIC partner potential to a repulsive Coulomb potential 
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for A = l,k = l,Z = 6, while Fig. 3b shows the corresponding wave functions. For 
comparison the original Coulomb potentials and wave functions are also shown dot- 
ted. It is seen that the potential which holds a bound state of positive energy shows 
an oscillatory behavior about the Coulomb potential, Vc, as is also evident from the 
form of Eq. for V. Since the oscillating component vanishes whenever uq vanishes, 
we have V = V a.t each node of Uq. Compared to the original, unnormalizable wave 
function, the BIC wave function in both cases shows a damped behavior due to the 
denominator function. This is also seen in the figures. 

A similar behavior is also expected for other radially symmetric potentials with 
a continuous spectrum of positive eigenvalues. For one-dimensional potentials, the 
situation is not so clear cut. Our method works for the Morse potential which is 
steeply rising on the negative x-axis with correspondingly damped wave functions. It 
also works for the case of a particle in a constant electric field for similar reasons. 
For potentials, such as V(x)=— Vq sech'^{x), the integral Iq Eq.(|^) is not convergent if 
the starting point is chosen at — oo, and it gets negative contributions if the starting 
point is selected at finite x-values. This leads to a vanishing denominator function in 
the expressions for some wave functions which makes them unacceptable. 

III. The Two Parameter Family of Potentials 

In Section II, we have seen that a straightforward procedure exists using the SUSY 
technique for generating a completely isospectral one-parameter family of potentials 
and that these potentials have a bound state in the continuum if we select as a starting 
point a positive energy solution of the Schrodinger equation for any potential V{t). 
We now show how this procedure can be extended to construct two-parameter fami- 
lies which contain two BICs. 

In constructing the new wave functions for the one-parameter family, Eq. (|^), we 
observe that the denominator function given in Eq. (Q) was all that was needed to 
create the BIC, while the operation in Eq. (|^) ensured that the wave functions for all 
the other states, there represented by ui, are a solution to the new potential. Note 
again, there is nothing special about the ordering of the two energy values nor the 
relative magnitude of Eq and Ei, therefore we can repeat this procedure by applying 
the theorem to the wave functions and the potential of the one-parameter family, but 
this time we transform the state at Ei into a BIC. The state at Eq, which already is 
a BIC, is transformed in the step of Eq. (j^), suitably modified, to become a solution 
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to the new potential. In this way we obtain the two parameter family of potentials 



V{r- A, Ai) = - 2[ln{h + Ai)]" = V 



h + Ai (Ji + Ai 



with the solutions of the corresponding Schrodinger equation 

Uq = {Eq - Ei)uo + ill W{ui, Mo), 
1 . 



and 



Ui 



h 



/i + Ai 



-til, 



(14) 

(15) 
(16) 

(17) 



The precise relationship of the new potential and its wave functions, which are now 
both BICs, is illustrated in the last column of the table. 



While the compact form of Eqs. ([T^ - |T^) explicitly shows the method of con- 
struction, it is useful to observe that the integral Ji can be conveniently re-cast into 
a simpler form which contains integrals of the form 

ul{r')dr\ (18) 
involving the original wave functions only. Making use of Eq. (||) for mi, we get 



h 



{Er-Eoful + -^^,+2{E,~E,) 



(/o + A) 



W 



dr'. 



The second term is integrated by parts as 



Un 



(/o + A)' 



W\r')dr' 



/o + A 



+ 



2WW' 
(Jo + A) 



dr'. 



(19) 



(20) 



We now use Eq. (^ for the derivative of a Wronskian of two solutions of the Schrodinger 
equation to rewrite the second term and observe, that it exactly cancels the last term 
in Eq. (|19D. We therefore have 



/o + A 



+ (^i-Eo)'/i(r). 



(21) 



Here we have made use of the fact that our boundary conditions imply that W(0)=0. 



As an example, we evaluate the two-parameter potential 



V = V -2 



ln{{Io + X)[{E,-Eofli- 



/o + A 



+ Ai] 



(22) 
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The argument of the logarithm can be rewritten as 



{El - EoYloh - W\r) + AAi + A(Ei - EoYh + XJo- (23) 

We happen to have transformed first the state at energy Eq into a BIG and then, 
in the second step, the state at Ei, which introduced the parameters A and Ai. Let us 
now consider applying our procedure in the reverse order, that is let us first transform 
the state at energy Ei into a BIG and then the state at energy Eq, producing the 
parameters /i and /ii. For this situation, the argument corresponding to Eq. p^) is 



{E, - E^fhh - W\r) + ^^^^^ + ^^^{E^ - E^fh + ^^h. (24) 
Glearly, one expects symmetry. This is guaranteed if the parameters are related by 

^i = A(^i - EqY (25) 

/ii = Ai/(^i - ^o)'- (26) 

This also leads to the same two-parameter wave functions. We also note that trans- 
forming any state twice by Eq. (§) does not create a second denominator or anything 
else new, but simply changes the value of the parameter A as shown in ref . []12| . Finally, 



relation (|T5D ensures that all other eigenstates will be solutions to the new potentials. 

Figures 4 and 5 illustrate the appearance of potentials and wave functions for 
representative choices of A and Ai. Glearly, various choices of the parameters A and 
Ai lead to quite different looking potentials V . We note in Fig. 5 that, as Ai grows 
large, the two-parameter potential approaches the shape of the one-parameter family 
shown in Fig. 1 and as discussed analytically above. For easy comparison with other 
works, we have chosen k\ = 2k in the figures, however, the ratio of Ei to Eq need not 
be integral. 



IV. Summary 



We have demonstrated how the SUSY method, originally conceived for discrete spec- 
tra, can be generalised for the construction of BlGs. We were able to show how to 
generate a one-parameter family of potentials which possess a localized positive energy 
state, starting from a potential V(r) which has a continuum of positive energy states. 
The only requirement that V(r) must satisfy in order that it have such a continuum is 
that it approaches a constant as r — oo. Then the solution of the Schrodinger equa- 
tion with the potential V(r) is, oscillatory at large r, which we can take to be of the 
form sin(kr). Therefore the integral Jq Eq.(^, will be of the form of a constant plus 
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J^^ sin'^{kr') dr', where tq can always be found such that, for r > ro the solution of 
the Schrodinger equation is approximately proportional to sin(kr). This means that 
-^0 = Ci + r /2 + sin{kr) /{4k), where Ci is a constant. Therefore, uo, Eq-@ will vanish 
at large r as 1/r, making it a normalizable state. Thus our procedure for constructing 
a BIG from an initial potential V(r) is actually valid for any spherically symmetric 
potential which approaches a constant as r oo. The situation is more complex 
for one- dimensional potentials as discussed in the text. The SUSY procedure has in 
common with the original Von Neumann- Wigner [Q] method that it makes the wave 
functions normalizable by generating a denominator function which grows with r as 
r ^ oo. In the case of V=0, our denominator function, containing Iq is a special case 
of the form used by Von Neumann and Wigner. We illustrated the one-parameter 
method for two interesting and analytically solvable cases: V(r)= 0, the free particle, 
and V(r)=Z/r, the Coulomb potential. The procedure was readily extended to obtain 
two-parameter families with two BICs at arbitrarily selected energies. 

It is pleasure to thank Prof. W. Y. Keung for many helpful discussions and for 
first bringing the existence of bound states in the continuum to our attention. This 
work was supported in part by the U. S. Department of Energy under grant DE- 
FG02-84ER40173. 
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Figure Captions 



Fig. 1 shows two examples of potentials V{r) (solid) and the associated BIC wave 
functions Uo{t) (dashed) in the one-parameter family starting from V(r)=0 for k=1.0. 
Fig. la is for small lambda (A — 0.5) and Fig. lb for large lambda (A — 5.0). 

Fig. 2a shows the BIC potential (solid) derived from the attractive Coulomb po- 
tential which is also shown for comparison (dotted). Observe how the BIC potential 
oscillates around the original Coulomb potential. 

Fig. 2b shows the corresponding BIC wave function (solid) and, for comparison, 
the original Coulomb wave function (dotted). The damping of the BIC wave function, 
which makes it normalizable is evident. 

Fig. 3a shows the BIC potential (solid) derived from the repulsive Coulomb poten- 
tial which is also shown for comparison (dotted). Again, the BIC potential oscillates 
around the original Coulomb potential. 

Fig. 3b shows the corresponding BIC wave function (solid) and, for comparison, 
the original Coulomb wave function (dotted). The damping of the BIC wave function, 
which makes it normalizable is evident. 

Fig. 4a shows V, a typical member of the two-parameter family of BIC potentials 
obtained from V(r)=0 for k=1.0, ki=2.0 and A = 1.0, Ai = 2.0. This potential sup- 
ports two bound states in the continuum at = 1, and at Ei — 4. The associated 
wave functions are shown in Fig. 4b, the lower state at Eq dashed, the higher one at 
El dotted. 

Fig. 5a shows V, a typical member of the two-parameter family of BIC potentials 
obtained from V(r)=0 for k=1.0, ki=2.0 and A = 1.0, Ai = 50. This potential also 
supports two bound states in the continuum at Eq = 1, and at Ei = 4. The associ- 
ated wave functions are shown in Fig. 5b, the lower state at Eq dashed, the higher 
one at Ei dotted. Because in this case Ai is relatively large, the potential approaches 
the shape of the one-parameter potential shown in Fig. la. Note that both energies 
Eq and El are above the maximum value of V. 
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Table. 



Potentials 



V = V - 2[ln(Io + A)]" 



V = V - 2[ln(ii + Ai)]" 



Wave functions 



El 



Ui 



Ui = (El - Eo)ui + UoW(uo, ui^ 



Ui 



ii+Ai 



Ui 



En 



Uo 



Uo = (Eo - Ei)uo + xiiW(ui, uq) 



Table Caption. 



The one-parameter family of potentials V{r; X) (central column) depends on the 
parameter A and has one bound state in the continuum at energy Eq with wave func- 
tion uq. Note that all other new states represented by Ui at Ei arc not normalizable. 
The right column shows the two-parameter family of potentials V{r; A, Ai), depending 
on the parameters A and Ai, which now has two normalizable states Uq and Ui in the 
continuum. Both families of potentials are generated from the non-normalizablc scat- 
tering states Uq and Ui of the original potential V(r) shown in the first column. Using 
the theorem described in the text, in the first step one produces a BIC at energy Eq 
and in the second step a BIC at energy Ei. While it is customary to denote the lower 
energy state by Eq, this is not necessary for our approach; Eq can also be above Ei. 
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